Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN I^TeX style file v2.2) 



Collapse and Fragmentation of Rotating Magnetized 
Clouds. II. 

Binary Formation and Fragmentation of First Cores 

! Masahiro N. Machida^* ,Tomoaki Matsumoto^f ,Tomoyuki Hanawa^J and 
O Kohji Tomisaka^§ 

^ Center for Frontier Science, Chiba University, Yayoicho 1-33, Inageku, Chiba 263-8522, Japan 
^ Faculty of Humanity and Environment, Hosei University, Fujimi, Chiyoda-ku, Tokyo 102-8160, Japan 
^ National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 



1—5 

0^ 



2 February 2008 



^ . ABSTRACT 

' Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized 

■<!;;j- I cloud are studied with use of three-dimensional MHD nested-grid simulations. After 

■ the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at 

' the center of the cloud. When the isothermal gas is stable for fragmentation in a con- 

, tracting disk, the adiabatic core often breaks into several fragments. Conditions for 

1^ ' fragmentation and binary formation are studied. All the cores which show fragmenta- 

I tions are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two 

r"| . patterns of fragmentation are found. (1) When a thin disk is supported by centrifugal 

^ i' force, the disk fragments through a ring configuration (ring fragmentation). This is 

Q , realized in a fast rotating adiabatic core as f2 > 0.2Tg^, where fl and rg represent 

^ • the angular rotation speed and the free-fall time of the core, respectively. (2) On the 

I other hand, the disk is deformed to an elongated bar in the isothermal stage for a 

Ci • strongly magnetized or rapidly rotating cloud. The bar breaks into 2-4 fragments 

ILJ ' (bar fragmentation) . Even if a disk is thin, the disk dominated by the magnetic force 

[ or thermal pressure is stable and forms a single compact body. In either ring or bar 

. fragmentation mode, the fragments contract and a pair of outflows are ejected from 

' the vicinities of the compact cores. The orbital angular momentum is larger than the 

5^ , spin angular momentum in the ring fragmentation. On the other hand, fragments of- 

ten quickly merge in the bar fragmentation, since the orbital angular momentum is 
smaller than the spin angular momentum in this case. Comparison with observations 
is also shown. 

Key words: binaries: general — ISM: jets and outflows — ISM: magnetic fields 
— MHD — stars: formation. 



1 INTRODUCTION 

Many stars are obser ved as m e mbers of binary or 
multi ple systems (e.g. lAbti Il983t IDuauennov fc Mavorl 
I1991I) . The bin ary frequency increa s es more in star- 
forming regions JCohen fc Kuhil Il979l: iDvck et al. 1 1 19821: 
IJov fc van Biesbroecklll994HRichichi et al. Il994l) . However, 
as summarized in IClarke fc Pringld Jl99jr the dynamical 
evolution of a small number of point-mass stars leads to a 
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system composed of one central binary or triple and many 
single escapers, which is a natural o utcome of a three- 
body encounter to make a close binary iBinnev fc Tremaind 
Il987l) . This does not explain the fact that binaries and 
multiples are common. Since the encounters in a cluster 
of stars or protostars cannot explain all the binary stars 
(Kroupa & Burkcrt 2001), we should devote attention to 
the direct binary formation through the fragmentation pro- 
cess i n the course of star formation (e.g. iBodenheimer et"ail 
I2OOOI). 



It is known that the molecular cloud collapses isother- 
mally until the gas density reaches n ~ 5 x 10^° cm"'' 
jLarsonI iTohlind 11982: iMasunaga fc Inuts uka 200(|. 

After the gas density exceeds this, gas becomes adiabatic. 
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because the cloud becomes optically thi ck against the dust 
thermal emissions and forms a first core jLarso^lIisi). The 
continuous accretion increases the mass of the first core and 
thus the gas density increases steadily in the core. Simul- 
taneously, the temperature rises steadily. Finally, molecu- 
lar hydrogen is dissociated after reaching T ~ lO'^K and 
n ~ 10^^ cm~^. Then the cloud begins to contract again as 
a result of the endoergic reaction. 

The fragmentation of molecular clo ud has been inves- 
tigat ed by many authors (for review, see iBodenheimer et all 
l2nnnr . Majority of researchers have focused on the isother- 
mal clouds. The initial conditions of the spherical isother- 
mal cloud with uniform density and rigid-body rotation are 
characterized by two parameters: thermal-to-gravitational 
energy ratio ao and rotational-to-gravitatio nal energy ratio 
00. T he evolution is divided into three types iMivama et al. I 
Il984 : Tsuribc & Inutsuka 1999). (1) a cloud with ao > I 
does not contract and oscillates. (2) a cloud with 0.5 < 
Qo ^ 1, which means the initial state is near the hydrostatic 
equilibrium from the Virial Theorem, experiences run-away 
collapse. The central density increases substantially in a fi- 
nite timescale. In this case, the contracting cloud does not 
fragment. (3) for a "cold" cloud such as Qq < 0.5, the cloud 
collapses in the direction of the angular momentum vec- 
tor and for ms a disk . The d isk is subject to fragme ntation 
as s hown in BonncU' ('1994), Bonnoll & Bate (1994a"b^ and 
IWhi tworth (1995) Although Cha & Whit worth (2003) have 
reported that even the cloud with ao — 0.6 fragments in the 
isothermal regime if it has a sufficiently strong differential- 
rotation, such a "warm" cloud does not fragment unless a 
ring is formed in the isothermal phase. 

Since a (nearly) rotation-supported disk is formed 
in the adiabatic phase, the disk fragments if it is suf- 
ficiently thin. Therefore, even if the gas does not frag- 
ment in the isothermal regime, it breaks into fragments 
after the adiabatic core has f ormed (the first core). Actu- 
allv. | Cha fc WhitworthI J2003ft . and iMatsumoto fc Hanawal 
i200a) assumed the barotropic equation of state for a gas 
and studied the evolution of a "warm" cloud with qq = 
0.6 — 0.76. They found that the thin disk fragments in the 
adiabatic phase even if ao is chosen for the cloud to expe - 
rience the run-away collapse. IMatsumoto fc Hanawal (l2003l) 
showed that a cloud with Qots ^ 0.05 fragments, where 
Qo and ta denote the initial central angular velocity and 
the free-fall time at the center (371/326*^0)^^^, respectively. 
By comparing different initial rotation laws they found that 
the central angular velocity is essential to whether the cloud 
fragments or not. This shows us that the fragmentation con- 
ditions are directly related to flottt rather than /3o. The 
thickness of the disk is also impor tant in the fragmenta- 
tion p rocess in the adiabatic stage (IMa tsumot o fc Hanawal 
l2003l: iMachida. Matsumoto fc Tomisakal |2004^ . Since the 
Toomore's Q value of the disk is smaller than 1 and the 
gravitationally most unstable wave-length is much smaller 
than the disk rad ius, the disk which fragment s is gravita- 
tionally unstable (IMatsumoto fc Hanawall2003ll . 

Disk structure is believed to be formed either by the 
rotation or the magnetic field. Almost all previous studies 
focus on the hydrodynamical contraction process. However, 
the molecular cloud is magnetized, and magnetic field lines 
alon g the major axis of molecular clouds are often obse rved 
fe.g. iTamura et al. lll995l : | Ward- Thompson et al. 1 1200(1 . In 



star-forming reg ions, the molecular outflows are seen to b e 
ubiquitou s (e.g. lOhashi et al |[l993 iBelloche et 1111200211 . 
iTomisakal (Il99d l200l |2002^ showed that these outflows 
are driven by the magnetic force in his two-dimensional 
axisymmetric magnet o hydro dynamical (MHD) simulations. 
iBasu fc Mouschoviaj iigg^T) pointed out that magnetic 
braking also plays an important role in a magnetically 
supported rotating cloud. Consequently, the magnetic field 
should play an important role in the star-formation process. 

When the cloud has no rotation motion, the magnetic 
field must have a positive effect to form a thin disk (see, 
e.g., Nakamura, Hanawa & Naka no|ll995l : [Tomisaka 199(|). 
However, in the rotating cloud, since the angular momen- 
tum can be removed from the core by magnetic force 
jBasu fc MouschoviasllTggi : ITomisakal I2OO0I) . the magnetic 
field may have a negative effect for fragmentation. To inves- 
tigate the effect of the magnetic field on fragmentation or 
binary formation, three-dimensional calculation taking ac- 
count of the magnetic field is needed. 

However, there have been few studies of these calcu- 
lations. It has not been cleared yet whether t he m a gnetic 
field promotes fragmentation or suppresses it. iBosd (l2002l) 
studied the evolution of magnetized molecular clouds and 
he claimed that the magnetic field promotes fragmentation. 
However, in his study, the angular momentum transfer by 
magnetic tension force is not included, because the mag- 
netic field is accounted for by an approximate form. Re- 
cently, iHosking fc WhitworthI (jgOO^) studied fragmentation 
of a magnetized cloud using their MHD SPH code. They 
found that fragmentation is prevented by the magnetic field 
because magnetic braking is so effective in their cloud. How- 
ever, they calculated only 4 models with simple initial con- 
ditions. See also Boss (2004) for a counterargument against 
iHosking fc WhitworthI (|2004^ . Thus, it has not become clear 
whether the molecular clouds observed can actually frag- 
ment or not. In our study, having a large parameter range of 
initial rotation speed and magnetic field strength, it is found 
that frag mentation is suppressed by th e magnetic field, as 
denoted in lHosking fc WhitworthI i2004ft . although the mag- 
netic braking is not so effective. 

Ambipolar diffusion affects protostellar collapse espe- 
cially at high density exceeding n ~ 10^^ cm~^. However, the 
diffusion t imescale has uncerta inty depending on the ioniza- 
tion rate iNakano et al. ' 2002). Furthermore, a calculation 
of ambipolar diffusion requires large computational costs: 
Hosking fc Whitworth (2004) calculate only a few models, 
and Boss (2002) adopts a simple approximation in the treat- 
ment of the diffusion. We therefore adopt the ideal MHD 
in order to perform an extensive parameter survey. Nakano 
et al. (2002) shows that the magnetic field is coupled with 
gas in n < lO^^"^'^ cm~^, indicating that the assumption 
of an ideal MHD is valid in the isothermal phase. In the 
adiabatic phase, the number density in the adiabatic core 
exceeds ~ lO^'^cm"^, and the magnetic field begins to de- 
couple. Our simulation may therefore overestimate the an- 
gular momentum transfer by magnetic field, especially in the 
dense fragments. 

We study the evolution and fragmentation of molecu- 
lar clouds using full three-dimensional MHD simulations. In 
this study, we use a nested grid code, which always main- 
tains sufficient spatial resolution in the central region. In 
this simulation, structures in the range of 5 orders of mag- 
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nitude in spatial extent (e.g. ~ 10® AU — lOAU) are resolved, 
which corresponds to 14 orders of magnitude in the density 
contrast (e.g. 5 x 10^ cm"^ - lO^^'cm"^). 

We calculated 144 models with different magnetic 
field strengths, rotation speeds, and initial amplitudes 
of non-axisymmetric perturbations. We have already re- 
ported a part of the result and a condition of fragmenta- 
tion in Machida, Tomisaka & Matsumoto (2004; hereafter 
MTM04). MTM04 shows that fragmentation occurs if the 
first core has a large oblateness (radius-height ratio). If the 
disk is thin enough (oblateness > 4), a non-axisymmetric 
mode develops in it after the gas becomes adiabatic. When 
a thick disk or a core is left at the core formation epoch 
(oblateness < 4), however, an axisymmetric core continues 
to contract. This seems to indicate single-star formation. If 
the non-axisymmetric perturbation has grown significantly 
and a bar appears in the isothermal stage, such a bar frag- 
ments in the following adiabatic accretion stage. 

Subsequent to the companion paper (Machida et al. 
2004; hereafter Paper I), we report results of three- 
dimensional MHD simulations in detail and give a condi- 
tion for binary formation. Paper I has studied the early evo- 
lution in which gas behaves isothermally. We have found 
that (1) a disk is formed in a direction perpendicular to the 
magnetic field and the rotation axis, and (2) the ratio of 
magnetic flux density to the square root of the gas pressure 
[B /{8n c'^p)^^'^] and that of the angular rotation speed to 
the gravitational free-fall rate [Q/{A-KGpy^^], both of which 
are evaluated at the center of the cloud, are well correlated 
at the core-formation epoch as 
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We evaluated the amplitude of the nonaxisymmetricity 
with the axis ratio viewed from the top at the end of 
the isothermal phase. This controls the subsequent evolu- 
tion. We calculated cloud evolution in the density range 



from n = 5 X 10 cm 



to 



In this paper. 



we present the evolution of the adiabatic accretion stage 
(n > 5 X 10^''cm~'^), while that of the isothermal collapse 
stage (n ^ 5 X 10^'^ cm~^) was shown in Paper I. Model 
and numerical method are given in §2, and the results of 
our calculation are presented in §3. In §4, we discuss the 
fragmentation conditions and the angular momentum redis- 
tribution between orbital and spin angular momenta after 
the fragmentation. 



2 MODEL AND NUMERICAL METHOD 

The initial model and numerical method are the same as 
those of MTM04 and Paper I. Here we describe them 
briefly. We assume that gas obeys the ideal MHD equations 
(eqs.[I.l]-[I.4]) for simplicity, where equation [I.l] represents 
equation (1) of Paper I. Interstellar gas behaves isothermally 
for n < 5 X lO^'' cm^'^ and adiabatically wit h the specific 
heat ratio of 7 = 7/5 for n > 5 x 10^° c m~^ iLarsonillQGgl: 
lTohlinell983 : lMasunaga fc Inutsukal200dl . Thus, we assume 
a barotropic equation of state as 
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where Cs denotes the isothermal sound s peed and Pcri = 
1.9 X 10~" g cm~^ or rtcri = 5 x lO" c m~^ jBonneU fc Batj 
'l994d; iBate. BonneU fc Brominl l2002h . As in MTM04 and 
Paper I, a filamentary cloud of infinite length is assumed as 
the initial state, which is in a hydrostatic balance radially. 
To this hydrostatic distribution we added axisymmetric den- 
sity perturbation 5pz{z) as well as non-axisymmetric density 
5ptp{r, ip) and magnetic flux density perturbations 5B^{r, ip) 
as equations [1.7] and [1.9]. 

We chose to describe the amplitude of axisymmetric 
density perturbation are proportional to the unperturbed 
density with the amplitude fixed at 10%. The initial model 
is characterized by three nondimensional parameters: the 
magnetic-to-thermal pressure ratio. 



= -B|c,o/(4vrpc,oCs,o), 



(3) 



the angular rotation velocity normalized by the free-fall 
timescale. 



oj = VLc,q/ (47rGpc,o)^^^ , 



(4) 



(2) 



and the initial amplitude of the non-axisymmetric pertur- 
bation, A^. The former two specify the equilibrium model, 
while the latter does the perturbations. We made 144 mod- 
els by combining the above parameters and calculated them 
with the 3D nested grid code, as we did in MTM04 and 
Paper I. We assume that the cloud has a central density 
of no = 5 X 10^ cm~^ initially, except for some models in 
which the cloud has a central density of 5 x 10* cm^'' and 
5 X 10® cm"''. In the following sections, we mention the evo- 
lution of the clouds only with no = 5 x lO'^ cm"'' and discuss 
its dependency on the initial density in §4.1. 



3 RESULTS 

When the central density exceeds ric — 5x 10^° cm"'', the gas 
temperature begins to rise, which supports t he gas and form s 
a static core (the first core) as denoted in ILarsoril il96Sfl . 
When the first core rotates or has a magnetic field, various 
substructures (c.f. core, disk, ring, and bar) appear. That is, 
after the gas becomes adia batic, fragmentation may occ ur as 
shown in MTM04 (see also lMatsumoto fc Hanawai2003l for a 
hydrodynamical case). In this section, we discuss cloud evo- 
lution in the adiabatic accretion phase (nc > 5 x 10*" cm"''), 
devoting particular attention to the fragmentation patterns. 
As shown in MTM04, the adiabatic core shows three char- 
acteristically different types of evolution. In the case with a 
strong magnetic field, magnetic braking works and the an- 
gular momentum is efficiently removed from the core. This 
results in a compact adiabatic core. From a fast rotating 
cloud, a thin disk appears and this evolves into a ring. The 
ring breaks through a non-axisymmetric spiral mode and 
leads to binary fragments, which is called the 'ring fragmen- 
tation' mode. The last case, in which a large amplitude of 
the initial nonaxisymetric perturbation forms a bar, even- 
tually fragments into two or several pieces. The fragments 
contract and form compact cores. This is called the 'bar 
fragmentation' mode. 

We show these three models, the core, ring fragmenta- 
tion and bar fragmentation models in greater detail in the 
following sections. The model parameters and physical quan- 
tities are summarized in Tabled In this table, the models 
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named BL (core model), CS (ring fragmentation model), and 
DL (bar fragmentation model) are the same as those in Pa- 
per 1. The ratio of the total to the magnetically critical mass, 
M/MB,cri ( MB,cri = -Bz,co/27rG^'^^ ) , is also described in this 
table. All models are magnetically supercritical cloud at ini- 
tial, because M/MB,cri exceeds unity. The evolution of these 
clouds in the adiabatic accretion phase {ric ^ 5 x 10^" cm"'') 
is shown in this paper, while that in the isothermal collapse 
phase (ric < 5 x 10^° cm~^) is described in Paper I. 

3.1 Core Model 

First, we show the typical evolution of the 'core' model 
(model BL) in the adiabatic accretion phase. A cloud with a 
weak magnetic field {a < 0.1 ) and slow rotation {uj < 0.1) 
forms a thick disk or an oblate spheroid at the end of the 
isothermal collapse phase (ric ^ 5 x 10^" cm"''). In these 
models, the central density continues to increase. In the 
central region, a dense and compact core appears, but this 
evolves without any indications of fragmentation during the 
adiabatic accretion phase {ric ^ 5 x 10^" cm""'). Figs.Qand 
|5| show the evolution of model BL (a, lj, A^p) — (0.1, 0.01, 
0.2). We calculate this model from the maximum density 
of 5 X 10^ cm~^ to 5 X 10^® cm"''. However, in this section, 
we focus on the evolution only in the adiabatic accretion 
phase [see §4.2 of Paper I for the isothermal evolution]. At 
the beginning of the adiabatic accretion stage, the central 
oblateness and axis ratio reach, respectively, eob ~ 5.3 and 
Ear = 0.23. ^ In Figs. and |5| panels (a) through (d) are 
snapshots at tc — 46 yr (a), 182 yr (b), 217 yr (c) and 
227 yr (d) after the core formation epoch, respectively. We 
define the core formation epoch (tc = 0) as the epoch in 
which the central density reaches ric = 5 x 10^*^ cm"''. The 
thick lines in Fig.0 which correspond to the contour lines of 
n = 5 X lO^'' cm"'', indicate the adiabatic core. In the lower 
panels of Fig. |21 other thick lines are drawn to indicate the 
outflow region with Vz > 0. 

Fig. (a) shows that the adiabatic core (inner thick 
line) has an almost axisymmetric shape and the diameter of 
the core is about 45 AU on z=0 plane. Gas inflows radially 
and the inflow speed reaches Vr = —0.54 kms"^ on the z —0 
plane, while the gas has a low rotation speed of = 0.06 
kms"^. At the same time, the gas forms a thick disk with a 
height of 36 AU as shown in the upper panel of Fig. |5|(a). 
It shows that gas flows vertically outside the disk, while it 
flows radially inside the disk [the lower panel of Fig.j^Ja)]. 

The central density increases from 9.1 x 10^° cm"'' 
[Fig. 111(a)] to 1.4 X 10^^ cm"^ [Fig. □(b)]. The first core 
changed shape and become ellipsoid, as shown in Fig. 0(b). 
This shows that the non-axisymmetry grows after the gas 
becomes adiabatic. Fig. |21 (b) (x—O cut) shows that this 
ellipsoidal core thickens and the vertical-to-horizontal axis 
ratio of the adiabatic core increases from 0.2 [Fig. |5| (a)] to 
0.4 because the gas temperature and, accordingly, the sound 
speed increases in this phase. This panel shows two shocks. 
The outer shock front is seen at 2; ~ ±40AU, while the inner 
shock is seen in the region of \z\ < 40AU and \y\ < 15AU. 

^ Let the length of principal axes ai ^ a2 in the x — y plane and 
in the 2— direction. The oblateness (Cob) a^nd axis ratio (ear) 
are defined y'ai 02/03 and 01/02 — 1 [see §4.1 of Paper I ]. 



Fig. |5Jb) shows that the gas inside the inner shock has a 
higher infall velocity {vz — —1.0km s"^) than that outside 
the inner shock {vz — —0.2km s"^). 

Fig. Wic) and (d) show that the adiabatic core is com- 
posed of a small core and a disk. Figs, □c) and (d) indicate 
that the outer boundary of the adiabatic core appears al- 
most the same as that of Fig. ^b) [Panels (c) and (d) are 4 
and 8 times more magnified than that of panel (b), respec- 
tively]. At the stage of Fig.0[d), the central density reaches 
1.1 X 10^'' cm""' and a dense compact core with a steep den- 
sity gradient has formed at the center. A disk with density 
in the range from 2 x 10^^ cm"^ to 5 x 10^^ cm"'' surrounds 
the core. 

As shown bv lTomisakal (Il998l . l2002h . magnetic field lines 
are tightly twisted after the adiabatic core formation. Fig. |H] 
left panel shows the magnetic field lines (streamlines) at the 
same stage of Fig. 0(d). This shows the magnetic field lines 
are tightly twisted by the rotating adiabatic core (red iso- 
density surface). The magnetic field lines begin to twist be- 
cause the collapse timescale exceeds the rotational timescale 
in the adiabatic accretion phase. The twisted field transfers 
the angular momentum of the adiabatic core to the outflow, 
and the outflow emerges from the core and brings the an- 
gular momentum from the core. In this model, the angular 
momentum is extracted through the gravitational and pres- 
sure torque since the adiabatic core is appreciable elongated. 
Thus, it is difficult to evaluate the angular momentum ex- 
traction due to the outflow separately. However, we have 
conflrmed that the angular momentum is largely removed 
by the outflow in the model of which parameters are the 
same as model BL except for = 0.01. The angular mo- 
mentum extraction is mostly due to the outflow since the 
core is nearly axisymmetric in the model of A^p — 0.01. For 
this reason, the angular momentum of the adiabatic core 
continues to be removed by the magnetic tension, and thus 
collapse continues. 

In the lower panels of Fig.jSJc) and (d), a strong outflow 
is seen, which is driven by the spinning of the adiabatic core. 
The outflow reaches z ~ 10 AU with Wz,out — 2.86 kms"^ of 
the maximum outflow speed at the stage of Fig.|5|(c), and ex- 
tends further to z ~ 18AU with fz,out — 2.90 km s^^ at the 
stage of Fig. 121(d). Fig. 2(d) upper panel indicates a typical 
structure of the adiabatic core (thick line) from which out- 
flow is ejected. There is a horn which extends from the oblate 
spheroidal adiabatic core. This separates the inflowing and 
outflowing gases [lower panel of Fig. 2 (d)]. Inside the horn 
including the rotation axis, gas is outflowing with a speed 
up to Vz — 2.9km s"^ and the outside the horn near the disk 
mid-plane gas is inflowing radially with Vr — — 0.5km s"^. 

This conflguration of a compact core with a disk is seen 
also in other models forming a single core in the adiabatic 
phase. Some hydrodynamical models with no magnetic field 
a — {ui < 0.02) belong to this core-type solution. How- 
ever, this core -I- disk configuration does not appear in the 
models without magnetic field. In such models of a = 0, the 
adiabatic core is a simple thin disk, and this adiabatic disk 
gradually contracts. This seems to indicate that the com- 
pact core in the adiabatic outer core or disk is composed of 
the gas from which angular momentum is extracted by mag- 
netic braking. The magnetic braking supresses the bar mode 
instability which took place in Bate (1998) and Matsumoto 
& Hanawa (2003), who assumed no magnetic field. 
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3.2 Ring Fragmentation Model 

In this section, we show the second typical evolution, e.g. 
'ring fragmentation' model (model CS) in the adiabatic 
ac cretion phase. This type of fragmentation is als o see n 
in iNorman fc WilsonI J 1978ft : Flbhline U98g}; Bosa JlQSCll . 
When the bar mode does not evolve so much in the isother- 
mal collapse phase (see §4.3 of Paper I), the central gas forms 
an axisymmetric disk at the beginning of the adiabatic ac- 
cretion phase. In the later adiabatic accretion phase, the 
disk is deformed to a ring, which is defined in this paper, 
as a structure where the central density is lower than the 
surroundings. 

Figs. 01 and |S] are the same as Figs. Q and |21 but for the 

'ring fragmentation' model CS [{a, uj, A^) = {0.Q1, 0.5, 0.01)]. 
In consequence to §4.3 of Paper I for the isothermal evolu- 
tion, we focus on the evolution in the adiabatic accretion 
phase. In this model, the gas interior to the isodensity sur- 
face of 1/10 of the maximum density has an oblateness of 5.1 
and the axis ratio of 0.2 (the major-to-minor axis ratio on 
the x-y plane is equal to 1.2) at the core formation epoch. 
This means that a thin and almost axisymmetric disk is 
formed in the isothermal collapse phase. Panels (a) through 
(d) are snapshots at =473 yr, 879 yr, 1344 yr, and 2134 
yr from the core formation in Figs. |1| and |K| 

473 yrs after the core formation epoch, the density peak 
moves from the center to the periphery (r ~ 20 AU); that 
is, the adiabatic core changes shape and becomes a ring. At 
the same time, an m = 2 non-axisymmetric mode grows and 
the core begins to fragment. In Fig. 21 (a), we can see two 
fragments whose density peaks are located at {x, y) = (15 
AU, 20 AU) and (-15 AU, -20 AU). The axis ratio of the 
adiabatic core is ear ~ 0.5 (major-to-minor a:xis ratio ~ 1.5). 
In Fig. 0(a), we find two nested disks, in which the outer 
disk extends in the region — 50AU ^ z < 50 AU, while the 
inner adiabatic disk (inside a thick line) extends near —5 AU 
< 2 < 5 AU and -40 AU < a; < 40 AU. Both the inner and 
outer disks are bound by the shock. The outer shock, which 
is similar to that seen in the upper panel of Fig. |5Jb), has 
occurred in the isothermal collapse phase (see §4.3 of Paper 
I). Gas falls vertically and has a maximum infall velocity 
of 1.5 kms~^ just outside the outer shock front [z ~ ±50 
AU). Passing through the shock front, gas decelerates to 
Vz — — O.lkms"^ of the infall velocity. However, the infalling 
gas near the 2-axis accelerates to 1.2 km s~^, similar to the 
core model of Fig. |21 

The low-contrast spiral structure seen in Fig.|ljb) grows 
with time. Moreover, the adiabatic core indicated by the 
thick line itself evolves into two fragments. In this pro- 
cess, a pair of adiabatic cores form. The peak density of 
the fragments is as large as 2.7 x 10^^ cm~^, and the peaks 
are located at {x, y) ~ (5 AU, -40 AU) and (-5 AU, 
40 AU). Comparing the locations of the density peaks of 
Fig.|l](a) and (b), we can see the fragments revolve counter- 
clockwise around the center and simultaneously move out- 
ward. On the other hand, the central density decreases to 
ric — 2 X 10^'^ cm~"^, and in this range, gas behaves isother- 
mally and no longer adiabatically. Consequently, a ring-like 
or two-arm spiral structure is established by this stage. 

The fragments and ring structure at tc = 1344 yr 
are shown in Fig. 0] (c). The density peaks of the frag- 
ments are located at [x, y) ~ (50 AU, —50 AU) and (—50 



AU, 50 AU), and the maximum density increases up to 
~ 10^^ cm""^, while the gas density at the center decreases to 
~ 5 X 10^ cm~'^. Thus, there is a density contrast of 1000:1 
between the ring or the spiral fragments and the center. 
This panel shows four adiabatic cores. The lower density 
adiabatic cores (n ~ 5 x 10^° cm~^) near {x,y) ~ (20 AU, 80 
AU) and (-20 AU, -80 AU) are formed in the remnant of 
the ring and at the trailing end of the dense adiabatic cores. 
Fig. |S] (c) upper panel shows that the shock front expands 
further vertically. 

The fragments are seen near {x, y) ~ (65 AU, 5 AU) 
and (-65 AU, -5 AU) in Fig.0](d). At this stage, the ring 
structure of r ~ lOOAU — 170AU is moving outward. This 
expansion of the ring seems due to the fact that the accreting 
matter has a larger specific angular momentum than does 
the ring and thus the accretion increases the specific angu- 
lar momentum of the ring. The fragments continue to move 
outwardly and at the same time become more compact and 
dense, while the density of the ring gradually decreases. Be- 
tween Fig.|3](a) and (d), the fragments revolve ~ 180 degrees 
around their orbits, although they revolve ~ 720 degrees 
until this calculation ends. The fragments around the or- 
bits revolve rapidly around each other (u^j.orb = 1.6km s~^), 
while they rotate themselves slowly (u,p,spn ~ 0.5km s~^). 
It should be noticed that the angular rotation speeds of the 
orbital motion Slorb are 1.5 times higher than those of the 
spin motion fispn- 

Now let's move on to the relation between fragmenta- 
tion and the outflow. In the lower panels of Fig. |^ (a) - 
(d), thick lines indicate the isovelocity contours of Vz — 0, 
which mean the boundaries between the infall and outflow 
gases. Gas is outflowing inside the thick lines. These pan- 
els show that the outflow region expands and the outflow 
accelerates with time. In the lower panel of Fig. |S](a), we 
can see two outflow regions, one near the shock front at 
2: ~ 50 AU and the other near the adiabatic core {z > 10 
AU). The maximum outflow speed is at most Uz.out — 0.27 
kms^^ at this time. At the stage of Fig. |3 (b), the maxi- 
mum outflow speed reaches Vz,out — 0.58 kms~^, while the 
maximum infall speed reaches ~ 1.5 kms~^. As the outflow 
region expands, the inner accretion shock front moves in 
the z-direction from z ~ ±50AU [Fig. I^a)] to z ~ ±70AU 
[Fig. Elb)] . As the adiabatic cores (fragments) move radially 
outward, the outflow region expands radially, too. The lower 
panel of Fig. EIc) shows the outflow is further accelerated 
and its speed reaches Uz,out — 0.7 kms~^, while the infall 
velocity remains ~ 1.5 kms~^ at maximum. Driven by a 
strong outflow, the shock fronts move upward to z ~ ±100 
AU. 

As the ring expands, the outflow region also expands ra- 
dially. Fig.|5Jd) shows that two types of outflow exit. That is, 
the outer outflows are connected to the ring structure and 
the inner ones are connected to the fragments. The outer 
outflows have a maximum speed of 0.84 kms~^ near {x,z) ~ 
(125 AU, 90 AU) and (-125 AU, 90 AU). The inner outflows 
are seen just above the spinning fragments near (x, z) ~ (80 
AU, 50 AU), and (-80 AU, 50 AU) and have a complicated 
shape. In Fig.|S| we can see the relation of the magnetic field 
lines and the outflow regions more clearly. In this flgure, the 
magnetic fleld lines (left panel) and the distribution of out- 
flow regions (right panel) are shown for the I = 12 subgrid 
(the box size is equal to 180 AU) from a bird's eye view 
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for the same epoch of Fig. 2] (d). This shows that almost 
all the magnetic field lines near the adiabatic core {r <68 
AU) are launched from the two fragments. This is natural, 
since at least 75 % of the mass in the adiabatic core is con- 
centrated in these fragments. The left panel shows that the 
magnetic field lines (streamlines) are tangled in the region 
above the remnant of the ring, which is indicated by density 
contours on z — plane, and are strongly twisted by the 
orbital azimuthal motion of the ring. Fig. |S] also shows that 
the two types of outflow are driven by the rotational mo- 
tion of the ring and the spin of the fragments, respectively. 
The outer outflow has a cylindrical shape and surrounds the 
smaller-scale outflows. The outflows connected to the frag- 
ments occupy the region interior to the ring-connected outer 
outflow, and are accelerated with core contraction similar to 
the outflow appearing in the core model (Fig. |^. However, 
the maximum speed of the outflow is slower than that at- 
tained in the bar fragmentation case. This may be related 
to the fact the angular rotation speed of the spinning of the 
fragments, fispin, is at most comparable to that of the or- 
bital motion, Sloibit, in the ring. Other models which show 
the ring fragmentation indicate similar evolutions: fragments 
move outwardly and two types of outflow co-exist. 

3.3 Bar Fragmentation Model 

In this section, we show the third type: the 'bar fragmen- 
tation' model (model DL). The cloud with a large initial 
amplitude of the non-axisymmetric perturbation and either 
a high rotation speed or strong magnetic fleld has changed 
shape from a disk to a bar in the isothermal collapse phase 
[see Paper I §4.4]. The bar grows from a disk and flnally 
fragments, as shown in Nakamura & Li (2002, 2003). Figs.Q 
and|H|are the same as Figs.Qand|21but for model DL [{a, uj, 
A^) = (1.0, 0.5, 0.2)]. See § 4.4 of Paper I for the isothermal 
collapse phase of model DL. This cloud has an oblateness 
of 5.2 and an axis ratio of 15.1 at the end of the isothermal 
collapse phase. A very thin bar has already formed by the 
core formation epoch. Four panels of Figs. Q and |5] corre- 
spond to snapshots at tc = 773 yr (a), 1033 yr (b), 1133 
yr (c) and 1235 yr (d). In model DL, fragmentation occurs 
immediately after the core formation. 

In Fig.|7](a), two density enhancements are seen in the 
bar near {x, y) ~ (20 AU, 35 AU) and (-20 AU, -35 AU). 
This panel shows that gas flows into the bar mainly perpen- 
dicular to the major axis of the bar on the 2 = plane, while 
inside the bar, gas flows to the center along the major axis 
of the bar. Fig. |5] indicates a cut by the a; = plane which 
is almost perpendicular to the major axis of the bar. Fig.|H| 
(a) shows that upwardly moving gas is already distributed 
above the bar. 

In Figs. |3b) andlHIb), the density of the fragment at 
{x, y) ~ (15 AU, 32 AU) and (-15 AU, -32 AU) is in- 
creased up to 5.9 X lO^'^ cm~^, and the fragments are grad- 
ually transformed into a sphere from an elongated bar. On 
the other hand, the central density of the bar ric decreases 
from 2 X lO" cm"^ [Fig.|ZI(a)] to 6 x cm"^ [Fig- 13 (b)]. 
This indicates that gas flows into fragments even from the 
center of the bar. The fragments, which are located near [x, 
y) ~ (16 AU, 30 AU) and (-16AU, -30AU) in Fig.E^c), are 
more dense and compact than those in Fig. (b) . At this 
stage, the central part of the bar is squeezed and narrowed. 



In Fig. |7| (d), the density at the center decreases to 
3 x 10^" cm""^, and the former bar-shaped adiabatic core is 
now separated into two parts, while the maximum density 
of the fragments reaches ~ 10^'' cm~^ . At this stage, the 
fragments are located at (a;, y) ~ (10 AU, 25 AU), and (—10 
AU, -25 AU), having moved from (16 AU, 30 AU), and (-16 
AU, -30 AU) of Fig. [7|(c). This means that the fragments 
are gradually approaching each other as they contract. 

The lower panels of Fig.|Hl[a) and (b) indicate that large- 
scale outflows are ejected both from the front and rear of the 
bar. This outflow has a velocity of 0.52 kms~^ at maximum, 
while the maximum infall velocity is equal to 1.41 kms~^. 
These flgures show that the bar is enclosed in this outflow 
gas and the top of the outflow reaches ~ 130 AU at this 
time. 

The lower panel of Fig.|2l[c) clearly shows that two types 
of outflows co-exist, that is, one is connected to the bar, 
and the other is connected to each fragment. The outer bar- 
connected outflows are already seen in Fig.|2l(a) and (b). On 
the other hand, the inner fragment-connected outflows seem 
to be driven by the spin of the fragment cores, and these 
outflows grow after the fragments contract sufficiently. The 
outer outflow reaches z ~ 200 AU, which has 0.8 km s~^ of 
the maximum outflow velocity, while the inner one reaches 
z ~63 AU with 3.48 kms"^ of the maximum velocity. It 
should be noted that the maximum outflow speed of the 
inner outflow is much higher than that of the outer. At this 
stage, since the infall velocity is as fast as 1.57 km s~^, 
the maximum outflow speed exceeds the infall one. Fig. |H| 
(c) upper panel shows that gas moves upwardly (v^ > 0) 
for r < 50 AU and a low-density gas (ric ~ 3 x 10* cm"'') 
begins to rise toward the adiabatic core. This is caused by 
the inner fast outflow, where gas infalls at tc = 1033 yr of 

Fig.ETb). 

A strong mass ejection is seen in Fig. |H| (d), which is 
related to the contracting core. This outflow is faster and 
denser than the outer outflow. This is due to the fact that 
the spin angular rotation speed of the fragments is much 
higher than that of the rotation of the bar. The top of the 
inner outflow reaches z ~ 120 AU, and the flow has a max- 
imum outflow velocity of 5.46 kms~^. The fragments have 
changed shape to compact spheres from elongated prolate 
spheroids. At the same time, the fragments contract and in- 
crease their spin speed. Due to this rapid rotation, magnetic 
fleld lines are tightly twisted and fast outflows are ejected. 
Since the two lobes of the fast outflow have the same rota- 
tional direction, further calculation shows that they seem to 
merge to form a composite outflow. On the other hand, the 
outer outflows survive without merging. 

Fig.|^shows the magnetic field lines (left panel) and the 
outflow regions (right panel). This flgure shows snapshots of 
the same epoch of Figs. |3d) and IsTd) . and the box scale is 
equal to 300 AU. In the left panel, we find that the mag- 
netic fleld lines in this region are anchored mainly to the 
fragments and they are tightly twisted by the spin of the 
fragments. When the bar is formed, the magnetic field lines 
run perpendicularly to the bar and are slightly twisted, as 
in Fig. 11 of Paper I. After the fragmentation, the magnetic 
fleld lines gather around the fragments as they contract. At 
the same time, the fragments increase their spin rotation 
speed and, as a result, the magnetic fleld lines are twisted in 
the region above the fragments. Conversely, tightly twisted 
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magnetic field lines are evidence of rapid rotation of the frag- 
ments. In the left panel of Fig. |^ we can see a bend in the 
magnetic field lines between the straight outer and twisted 
inner field lines. Gas near this bend coincides with the out- 
flow appearing in the right panel of Fig.|^ Near this bend, 
the outer outflows which appear in the lower panels of Fig.|Hl 
are ejected. These outflows seem to be driven by the rota- 
tion of the bar remnant. The inner two outflows connected 
to the fragments are seen in the right panels of Fig. |U] It 
should be noted that the outflow region coincides with the 
prominent toroidal magnetic fields. 

The bar rotates very slowly as seen in Fig. |7] (a) - 
(d), while the ring rotates rapidly in the ring fragmenta- 
tion model [Fig.|l|(a) - (d)]. Owing to this slow rotation or 
less efficient centrifugal force, the fragments approach each 
other gradually in the bar fragmentation case. Although we 
have not seen the fragments merge in this model, bar frag- 
mentation models often show mutual merges and form a 
compound core (see Table |2j. 



4 DISCUSSION 

4.1 Cloud Fragmentation 

In this section, we discuss whether the fragmentation oc- 
curs or not, depending on the initial magnetic field strength 
(a) and the angular rotation speed (uj). In order to inves- 
tigate the fragmentation condition of the molecular cloud, 
we have studied 144 models with different parameter sets. 
Fragmentation occurs in 76 out of the 144 models studied. 
Ring fragmentations are seen in 53 of the 76 fragmenta- 
tion models, while bar fragmentation is seen in 23 models. 
The model parameters and results of all models are sum- 
marized in Table 121 In this table, the oblateness, axis ratio, 
central magnetic field strength and angular rotation speed 
at ric = 5 X 10^° cm~'^ are shown. The fragmentation mode 
(the core model, the ring fragmentation, or the bar fragmen- 
tation) and the final fate of the fragments, i.e. mergers or 
survivals, are also described. 

4.1.1 Evolution of the Isothermal Collapse Phase and The 
Magnetic Flux and Spin (B-Q,) Relation 

Although the clouds fragment only in the adiabatic accre- 
tion phase (ric ^ 5 x 10^" cm~'^), the cloud evolutions in the 
isothermal collapse phase (wc < 5 x 10^" cm~'^) are impor- 
tant in understanding the following fragmentation process in 
the adiabatic phase. The evolutions of the molecular cloud 
with different initial magnetic field strengths and angular 
rotation speeds in the isothermal collapse phase are studied 
in Paper I. Here, we summarize these evolutions briefly. 

In Paper I, we found the magnetic flux density and the 
angular rotation speed at the cloud center converge to a re- 
lation of equation Q in the process of isothermal collapse, 
irrespective of the initial amplitude of the non-axisymmetric 
perturbations. The first term of equation Q indicates the 
ratio of rotational to gravitational energy, and the second 
term means the ratio of the magnetic to thermal energy 
in the vicinity of the cloud center. Equation Q is plot- 
ted with a thick line in Fig. 1101 We named this curve the 
magnetic flux-spin [B-Q.) relation. As discussed in Paper 



I, there are four different evolutional patterns (c.f. models 
A, B, C and D) according to the initial magnetic flux den- 
sity and the rotation speed. Four arrows (models A, B, C 
and D), which meet the _B-f2 relation at the endpoint in 
Fig. 1101 represent these evolutional paths from the initial 
state (no = 5 x 10^ cm""^) to the beginning of the adia- 
batic accretion phase (tIc = 5 x 10^" cm~^). The evolutions 
of models B (the core model), C (the ring fragmentation 
model) and D (the bar fragmentation model) for the adi- 
abatic stage have been described in the preceding section. 
Since model A [(a, uj) = (0.01, 0.01)] forms a compact and 
dense core and shows evolution similar to that of model B 
in the adiabatic stage, we abbreviate this evolution in this 
paper. Models A-D have different initial amplitudes of the 
non-axisymmetric perturbation (for detailed evolution, see 
§4 of Paper I) . 

Since the B-fi space is divided into two by the above- 
mentioned B-Q. relation curve, the evolution of the cloud 
is also divided into two: that is, in the models distributed 
in the area below the B-Q. relation curve, the points move 
towards the upper-right or right, which means the cloud 
collapses spherically (see §4.1 and 4.2 of Paper I). In this 
region, the forces supporting the cloud radially (centrifugal 
and magnetic force) are insufficient relative to the gravity 
and the pressure force. On the other hand, in the models 
above the B-Q relation curve, the points move towards the 
lower-left, which means the cloud collapses vertically mak- 
ing a disk (see §4.3 and 4.4 of Paper I). Thus, we call the 
region in the lower or left side of the B-Q relation curve the 
'spherical collapse region' and that to its upper or right side 
as the 'vertical collapse region'. The magnetic flux density 
and angular rotation speed at the beginning of the adiabatic 
accretion phase (ric = 5 x 10^" cm~'^) are concentrated near 
or on the B-Q relation curve. We call models concentrated 
near the horizontal part of the B-Q relation curve (0 < 
Szc/(87rc?Pc)^''^ < 0.3 and Q^/{'^-kGpcY'^ - 0-2) 'rotation 
dominated disks', because the disk is mainly supported by 
the rotation. On the other hand, models near the vertical 
part of the B-Q relation curve (B,c/(87rc2pc)'''^ > 0.3 and 
< nc/(47rG'pc)''" < 0.2) are named 'magnetic dominated 
disks', because the disk is supported mainly by the magnetic 
held. 

4.1.2 Fragmentation Condition 

Fig. 1101 shows how the cloud fragments during dynamical 
contraction. The types of fragmentation are indicated by 
large symbols, as □(ring fragmentation), A (bar fragmenta- 
tion), 0(riiig or bar fragmentation, depending on A,p and 
ricfi), and x(no fragmentation). The large symbols appear 
at positions which represent the initial conditions of Qc and 
Bzc. Models indicated with Q fragment through the bar 
mode when the models have large non-axisymmetric per- 
turbations A^, or low initial densities nc,o, while they frag- 
ment through the ring mode when A,p is small or nc,o is high. 
From the distribution of symbols we can divide the param- 
eter space into several regions. Colored regions indicate the 
regions where the cloud fragments in the adiabatic accretion 
phase. Their colors denote the modes. Blue, red and violet 
regions represent, respectively, ring fragmentation, bar frag- 
mentation and either bar or ring fragmentation, correspond- 
ing to the amplitude of A^. 
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Small symbols, o (fragmentation) and + (no fragmen- 
tation), are plotted at the positions which represent the 
central magnetic field strengths (-Bzc) and angular rotation 
speeds (flc) at the beginning of the adiabatic contraction 
phase (5 x fO^" cm"''). Open and filled o symbols, indicate 
types of fragmentation as ring (open) and bar (filled). One 
large symbol (Q, □ or x) represents several models with 
the same magnetic flux density and angular rotation speed, 
but a different amplitudes of non-axisymmetric perturba- 
tion [A^ = 10""^, O.Of , 0.2, and 0.3), and central density 
{nc,o = 5 X fO^ cm-^ 5 x fO" cm"^ and 5 x 10*^ cm"^). There 
are 63 large symbols in this figure, while there are 144 small 
symbols. 

The ring fragmentation occurs in a large area in the 
initial parameter space shaded blue, while the bar fragmen- 
tation occurs when the initial angular velocity is larger than 
0.2 [i.e. ll) > 0.2] irrespective of the initial magnetic flux 
density. This is due to the fact that the evolution of the 
cloud differs greatly between the lower and the upper side of 
the B-n relation curve. That is, the axis ratio grows slowly 
in models below the B-Q relation curve. As discussed in 
MTM04 and Paper I, the non-axisymmetry begins to grow 
after a disk is formed, namely, growth of the axis ratio de- 
pends on the degree of oblateness. The oblateness grows in 

1/2 

proportion to pc in the area above the _B-f2 relation curve, 
while it grows to pl^^ below the curve. After the oblateness 
reaches 1-2, the axis ratio grows in proportion to pl^^ as 
shown in Paper I. Therefore, if the initial model parame- 
ters are chosen in the lower region, the non-axisymmetry 
does not grow sufficiently in the isothermal collapse phase 
because the oblateness grows only slowly. Even in this re- 
gion, the non-axisymmetry would grow if the initial cloud 
density were substantially lower than that of the observed 
cloud (ric <C 5 x 10^ cm~^). This is because in such a cloud, 
the isothermal collapse phase would continue for a long time. 
Thus, a cloud with parameters below the B-Q relation curve 
does not form an elongated bar. Therefore, only ring frag- 
mentation occurs for these parameters. 

We calculated models with = 10"^, 0.01, 0.2, and 
0.3 of the initial amplitude of the non-axisymmetric pertur- 
bation. In models in the lower region, no elongated bar ap- 
pears, even for large non-axisymmetric perturbation A^p = 
0.3. For example, models AL and BL, which have large ini- 
tial non-axisymmetric perturbations (A^p = 0.2), only have 
small axis ratios 0.068 and 0.23 at the end of the isothermal 
collapse phase (iic = 5 x lO^'' cm~^). However, if much larger 
non-axisymmetric perturbation (e.g. A^p > 0.5) is added to 
the model in the lower region, a bar may form in the isother- 
mal collapse phase. On the other hand, in the model whose 
parameters exist in the upper side of the curve, an elon- 
gated bar appears in the isothermal collapse phase, and the 
axis ratio grows even in the isothermal phase. Therefore, 
these models result in either a ring or bar fragmentation, 
depending on the initial amplitude of the non-axisymmetric 
perturbations. An area in which only bar fragmentation is 
possible appears in the upper right corner of the figure. In 
this region, the ring fragmentation is not possible even if 
a disk is formed in the isothermal collapse phase. The disk 
deforms to a compact core, because the angular momentum 
of the disk is removed by the magnetic tension. Neither ring 
nor bar fragmentation occurs in the unshaded area, in which 



a small core is formed because the initial rotation speed is 
slow or magnetic braking is effective. 



4-1.3 Ring Fragmentation 

The small open diamond symbols (ring fragmentation), 
which denote the value of the magnetic field strength and 
angular velocity at ric = 5 x 10^" cm"'^ in Fig. 1101 are dis- 
tributed in a narrow band of 0.2 < Q.^/ [A-KGpcf^'^ < 0.3. 
This means that whether the core displays the ring frag- 
mentation or not is determined only by the angular rota- 
tion speed at the stage of ric — 5 x 10^" cm"'', irrespec- 
tive of the magnetic field strength (see also Fig. IllH . If 
the cloud evolves to have a sufficiently fast rotation speed 
[fic/(47rGpc)^^'^ ~ 0.2] in the isothermal collapse phase, the 
ring fragmentation occurs in the subsequent adiabatic ac- 
cretion phase. As a result, the ring fragmentation condition 
could be denoted as follows: 

t7c/(47rGpc)'/' > 0.2, (5) 

at the beginning of the adiabatic accretion phase (ric = 
5 X 10^° cm"''). This condit ion is equivalent to that of 
iMatsumoto fc Hanawal i2003l) . in which they did not in- 
clude the magnetic field. If this condition has been satis- 
fied in the isothermal collapse phase, the cloud can frag- 
ment through a ring shape as shown in Figs. |1] and |21 The 
cloud in the area below the _B-f2 relation curve evolves to 
the upper-right at approximately 45 degrees in Fig. 1101 and 
the cloud approaches the B-Q. relation curve, as denoted in 
Paper I. Clouds with parameters in the upper region con- 
tinue to evolve to the lower-left until they reach the line. 
Therefore only clouds with initial parameters in either blue 
or violet areas can arrive at i7c/(47rGpc)^^^ = 0.2 and frag- 
mentation occurs. In other words, models with large a meet 
the magnetic-dominated part of the B-Q. relation where 
nc/iA-KGp^f^ < 0.2. 

The oblateness (sob) and the angular velocity 
[f2c/(47rGpc)^^^] are plotted at the end of the isothermal 
phase in Fig. llll (a'l. In this panel, symbols o, o and -I- repre- 
sent the ring, the bar and no fragmentation. As for ring frag- 
mentation models (o), the ring fragmentation occurs only in 
the upper-right area of Eoh > 3 and r2c/(47rGpc)^''^ > 0.2 
in Fig. 111! fa). The ring fragmentation does not occur with- 
out a thin disk even if a cloud has f2c/(47rGpc)^^^ > 0.2 in 
the isothermal collapse phase. The models in the upper-left 
region [sob < 3 and nc/(47rGpc)^/^ > 0.2] of Fig.IHIa) rep- 
resent the clouds which rotate rapidly and have thick disks 
or spheres. These clouds are formed if the initial density 
is high, because for such clouds, the isothermal contraction 
phase ends before a sufficient thin disk is formed. The oblate- 
ness at the end of the isothermal phase depends on the ini- 
tial density, because the oblateness increases with the cloud 
density as shown in Fig. 3 of Paper I. Thus, the clouds hav- 
ing extremely high densities (tIc > lO" cm"^) initially are 
distributed in this region. 

On the other hand, the clouds in the lower-left re- 
gion [£ob < 3 and nc/(47rGpc)^''^ < 0.2] of Fig. [TTIa) do 
not fragment through the ring mode either, in which these 
clouds initially have strong magnetic fields. These clouds 
form magnetic-dominated disks. In the lower-left region, the 
core appears as a thick disk or a sphere rotating slowly and 
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does not fragment. As a result, only the rotation dominated- 
disks in the upper-right region fragment through the ring 
mode. 

Fragmentation depends on the initial density especially 
for models in the area inside the B-Q, relation curve. We have 
plotted the area of the ring fragmentation (i.e. blue shaded 
area) assuming the initial cloud density to be 5 x lO'^ cm~^. 
However, if we adopt a higher initial density such as nc,o ~ 
5 X 10*cm~^ or nc,o = 5 x 10^ cm~^, the lower limit of 
the angular speed for fragmentation is raised because a 
only short period of the isothermal collapse phase is avail- 
able. In Fig. IIUI three models with different initial densi- 
ties (nc,o = 5 X 10^ cm"^, 5 x 10* cm"^, and 5 x 10*^ cm"^) 
are plotted. Comparing models with nc,o = 5 x 10^ cm~"^ 
and Tico = 5 X 10*cm~'^, the lower limit of the initial ro- 
tation speed for the ring fragmentation increases slightly 
from uj > 0.02 for nc,o = 5 x 10^ cm"^ to lj > 0.04 for 
nc,o = 5 X 10* cm~^ for a = 0. Because the model with a 
higher initial density does not reach the B-Q relation curve, 
the fragmentation condition, equation (O, is not fulfilled 
in the isothermal collapse phase for higher initial densities. 
Thus, although the lower limit of the initial rotation speed 
for the ring fragmentation depends on the initial density, 
this effect seems small. On the other hand, for models in 
the upper side of the B-Q relation curve, their evolutions 
hardly depend on the initial density. These clouds evolve 
rapidly and reach the B-Q relation curve when the central 
density increases for two orders of magnitude. Thus, these 
clouds induce fragmentation except for extremely high den- 
sity clouds (for example Uc > 10* cm~'^). 

As for model with q = 0, fragmentation does not occur 
for Lu < 0.02. However, if the initial density is lower than 
5 X 10^ cm"'', or in other words it has a long isothermal 
collapse phase, fragmentation may occur even for a cloud 
with UJ < 0.02. On the other hand, for a > 0.01, clouds 
with UJ < 0.02 never fragment, irrespective of their initial 
densities. This seems due to the fact that these clouds must 
evolve into magnetic dominated disks. 

4- 1-4 Bar Fragmentation 

Bar fragmentation occurs in the cloud with a large mag- 
netic field and an initially fast rotation speed. Small filled 
diamonds indicate the cores of the bar fragmentation in 
Fig. 1101 These symbols are almost distributed in the area 
of f2c/(47rGpc)^''^ > 0.2 where the symbols of ring fragmen- 
tation (open small diamonds) are also distributed. These two 
fragmentation modes co-exist in the same region. However, 
several bar fragmentation models are distributed outside the 
ring fragmentation region. Bar fragmentation occurs even in 
the models of strongly magnetized clouds (a > 1) if clouds 
have an initial rotation of uj ^ 0.3, while ring fragmentation 
occurs only in the rapidly rotating and weakly magnetized 
cloud. In the bar fragmentation models, fragmentation oc- 
curs in the elongated bar formed in the isothermal collapse 
phase. Fig. Illf b'l shows that six bar fragmentation models 
with angular speeds lower than f2c/(47rGpc)^''^ < 0.2 have 
large axis ratios Ear = 3.1, 3.3, 6.0, 6.8, 7.2 and 7.5. They 
are special models, that is, models initially in the region of 
Qc,o/i4:TvGpc,o)^^^ > 0.2 and a ^ I. Since they evolve to- 
ward the lower-left in the isothermal phase, nc/(47rGpc)^''^ 
becomes smaller than 0.2. 



Bar fragmentation occurs in the cloud with a large axis 
ratio, irrespective of the magnetic field strength and the an- 
gular speed at the end of the isothermal phase. For the bar 
formation, the cloud needs to be located initially above the 
B-Q, relation curve (i.e. in the vertical collapse region). The 
evolution of the axis ratio depends on the initial magnetic 
field strength and the rotation speed, and the axis ratio ob- 
tained at the end of the isothermal phase is an increasing 
function of a and uj as shown in Paper I. However, the cloud 
which is magnetized strongly or rotating slowly tends to 
form a small core. Such a core does not deform to an elon- 
gated bar and does not induce fragmentation. This indicates 
that the bar fragmentation region in Fig. 1101 resembles the 
ring fragmentation region above the B-Q relation curve. In 
this calculation, we added 30% of the non-axisymmetric per- 
turbation at maximum to the initial cloud. If we add much 
larger non-axisymmetric perturbations, the lower boundary 
of the bar fragmentation area (blue area) may extend down- 
ward in Fig. 1101 

4- 1.5 Comparison with Previous Work 

In this section, we compare our fragmentation conditions 
with previous work. Firstly, we consider the fragmenta- 
tion conditions of non-magnetized cloud. In our study, 
the fragmentation condition is denoted as uj > 0.02 for 
Wc.o = 5 X 10^ cm""^ (see Fig . IIOL in a cloud with a = 0. 
iMatsumoto fc Hanawal i2003l) studied the evolution of the 
non-magnetized rotating cloud with n = 2.6 x 10* cm~^ of 
initial density, and found that fragmentation occurs in the 
cloud with UJ > 0.026, which is mostly in agreement with 
our conditions for the non-magnetized cloud. 

T he fragment ation of the magnetized clouds is stud- 
iedbyHosi ((iooi) and lHosking fc WhitworthI ll2004l l . iBos j 
<2002l) assumes the cloud have 200 /xG (a ~ 4) magnetic 
field strength and n = 6 x 10* cm~^ initial density, and he 
adopt the angular rotation speed in the range of 10~^* s~^ 
- 10"^''* s"^ {uj = 0.01 - 0.12). It should be noted that he 
assumes the magnetic field strength to decrease with time 
to take into account the effect of the ambipolar diffusion as 
B = Bo,i{l—t/TAD), where tad ~ Wtff. As a result, he finds 
that all of the clouds fragment and concludes that the mag- 
netic field promotes the cloud fragmentation. These clouds 
are located to the right hand of the B-Q relation curve in 
Fig. IIUI In our study, the clouds in this region cannot frag- 
ment, because these clouds form compact and dense core in 
the adiabatic stage owing to magnetic braking. In his study, 
since magnetic braking is not taken into account, the specific 
angular momentum of infalling gas seems relatively high. 
Thus, only the effect which promotes cloud fragmentation 
(magnetic pressure and dilution of gravity by magnetic ten- 
sion) is included in his calculation. This is consistent with 
the fact that fragmentation occurs only in the cloud with 
fast rotation. 

On the other hand, a magnetized (764^G; a ~ 25) and 
non-magnetized (O^G; a — 0) clouds with the same angular 
velocity {uj ~ 0.2) and the same density (n = 2.8 x 10*^ cm"^) 
are investigated in Hosking & Whitworth ( 2004) . Note that 
they included both the ambipolar diffusio n and magn etic 
braking. We did not include the former while IBos3 l|2003) did 
the former in the approximation form but not the latter. The 
non-magnetized cloud fragments in iHosking fc WhitworthI 
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i2004l) . while the magnetized cloud does not. Comparing 
their parameters with ours, their magnetized cloud (a = 25, 
uj = 0.2) is in the non-fragmentation region in Fig. 1101 while 
their non-magnetized cloud {a = 0,uj = 0.2) is in the frag- 
mentation region. Therefore their result is consistent with 
our fragmentation conditions. We conclude that the mag- 
netic field suppress es the cloud fragmentation. The opposite 
result obtained by iBosd (l2002l) is likely to be ascribed to 
missing magnetic braking. 

4.2 Angular Momentum Redistribution 

In this section, we discuss how the angular momentum is re- 
distributed between the orbital and spin angular momenta 
in the typical bar and ring fragmentation models. After frag- 
mentation, the angular momentum of the host cloud is di- 
vided into orbital and spin angular momenta of fragments. 
As shown in >li-{.2l and the fragments seem to have large 
orbital but small spin momenta in ring fragmentation, while 
they have small orbital and large spin angular momenta in 
bar fragmentation. To estimate the angular momenta of the 
fragments quantitatively, we calculate the total (jtot), or- 
bital (jorb) and spin (jspn) specific angular momenta as fol- 
lows: 

jtot = jorb+ispn, (6) 

jorb = (rf X Vf) ■ ez, (7) 
i=pn = [ P(r) {(r - rf) x [v(r) - Vf]} dV ■ e./M(8) 

P>P25 

M = f p{r)dV, (9) 

P>P25 

where rf = J ^ p{r)rdV/M (J represents that 
summation should be done for one fragment), vt = 

/p>p P^^^ '^i^) dV/M, p25 and Bz mean respectively the po- 
sition vector of the mass center of the fragment, the bulk 
velocity of the fragment, a typical density of the fragment 
as p25 = 2.5 X 10^^ cm~^ and a unit vector of the z direc- 
tion. We identify a fragment as a gas which has a density of 
n ^ 2.5 X 10^^ cm~^ in equations @, Q, and ©. Although 
we have compared how the values depend on the choice of 
this density, 5 X 10^^ cm"^, 2.5 x 10^^ cm"^ and 5 x 10^^ cm"^ 
give essentially the same result. 

Fig. ll2l shows the three specific angular momenta (total, 
orbital and spin) of the fragments plotted against the time 
after fragmentation (tf). Three ring- fragmentation models 
(upper panels) and three bar-fragmentation models (lower 
panels) are plotted in this figure. Magnetic field strengths 
increase from the left to the right as a=0.01 [left panels: 
(a) and (d)], 0.1 [middle panels: (b) and (e)] and 1.0 [right 
panels: (c) and (f)]. The distances between the fragments 
are also plotted in each panel. The contour plots are inset in 
each panel to show the density distribution of 2=0 at about 
1200 yr after the fragmentation. 

4-2.1 Ring Fragmentation 

The upper panels of Fig. 1121 represent the angular momen- 
tum distribution in the ring fragmentation models. From the 
insets, the ring structure is still seen outside the fragments. 



In panel (a), which corresponds to a model with a weak mag- 
netic field, the fragment has 2.2 x 10^*cm^s~^ of the spin 
specific angular momentum, while it has 4.1 x lO^^cm'^s"^ 
of the orbital one at 200 yr after the fragmentation. The 
orbital angular momentum is about 20 times larger than 
the spin angular momentum. Fragment mass grows up to 
5.5 X 10"^ Mq in 1400 yr. At the same time, although the 
spin angular momentum increases up to 8 x lO^^cm'^ s~^, it 
is still 10 times smaller than the orbital angular momentum. 
The distance between the fragments has increased from 55 
AU (t/ = 200 yr) to 160 AU {tf = 1400 yr), because of 
the large angular momentum of the infiow gas. A similar in- 
crease in the spin angular momentum is observed in panels 
(b) and (c). It should be noted, in any case of ring fragmen- 
tation, the system has a large orbital angular momentum 
and a small spin momentum. The ratio of spin to the or- 
bital angular momenta is about 0.1, and this ratio seems to 
somewhat depend on the magnetic field strength. We can 
see from these panels that the total angular momentum is 
mainly composed of the orbital angular momentum. Tightly 
twisted magnetic field lines surround the ring. The magnetic 
field lines anchored to the fragments are weakly twisted as 
shown in Fig. |S| In panels (b) and (c), the distances be- 
tween the fragments decrease or oscillate, while they begin 
to increase in a further calculation {tc > 1400yr). For all the 
ring fragmentation models including those in panels (a) , (b) 
and (c), fragments survive without merging. The distance 
between fragments decreases in panel (b) and (c), while it 
continues to increase or keep constant without decreasing in 
further calculations. Therefore, these models seem suitable 
to form a binary system. 

4-2-2 Bar Fragmentation 

The lower panels of Fig. ll2l indicate the angular momentum 
redistribution in the bar fragmentation models. The inset 
contour plots in panels (d), (e) and (f) show that the frag- 
mentation occurs in the bar in these models. In panel (d) 
which is the case of a weakly magnetized cloud, the fragment 
has l.lx lO^^cm^ s~^ of the spin angular momentum, while it 
has 4.1 X lO^^cm^ s~^ of the orbital one at 200 yrs after the 
fragmentation. In this model, the orbital angular momen- 
tum is larger than the spin momentum, similar to the ring 
fragmentation models. However, the ratio of orbital to spin 
angular momenta is not too large compared with that of the 
ring fragmentation. After the fragment mass grows to about 
0.02 Mq {tf = 1000 yrs), the orbital angular momentum os- 
cillates while the spin momentum continues to increase. At 
the same time, the distance between the fragments decreases 
gradually. In this model, two fragments merge at about 10^ 
yrs after the fragmentation epoch. In the models shown in 
panel (e), the spin angular momentum continues to increase 
and catches up with the orbital one at 1260 yrs, and then 
orbital angular momentum begins to decrease rapidly, while 
the spin momentum decreases gradually. In a model with 
strong magnetic field [panel (f)], the spin angular momen- 
tum also catches up with the orbital one at 540 yrs, and then, 
both momenta increase. For models shown in panels (d) , (e) 
and (f), the spin angular momentum exceeds or is equal to 
the orbital one, when the fragment mass becomes sufficiently 
larger compared with that of the fragmentation epoch. Ow- 
ing to the decrease of the orbital angular momentum, the ra- 
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dial distance between the fragments decreases. Although the 
fragments do not merge in the models shown in (e) and (f ) in 
our calculation, the radial distance decreases steadily. Thus, 
the fragments seem to merge, if we continue the calculation 
further. Bar fragmentation models have the common feature 
of the angular momentum being equivalently redistributed 
into the orbital and spin angular momentum. Although the 
bar fragmentation occurs in 23 models in Fig. 1101 the mutual 
mergers between two fragments occur in 9 models. Ifowever, 
this number must increase for a longer simulation. On the 
other hand, the strong outflows are seen near the fragments 
in the bar fragmentation models as shown in Fig.|51 because 
of their large spin angular momenta. 

In conclusion, we show that the angular momentum re- 
distribution depends strongly on the fragmentation type, 
ring and bar, although it depends weakly on the magnetic 
field strength compared with panel (a)-(c) or (d)-(f). 

4.3 Comparison with Observation 

In this section, we apply our fragmentation conditions to 
several typical molecular clouds in which physical parame- 
ters are relatively well fixed. We find that the fragmentation 
region m Fig. ESI is weU fitted by an expression using initial 
magnetic field strength (B) and angular rotation speed (0.) 
as 

/ \ 1/2 

n B \ 



-f0.02 log 



2 X 10-22 



gem 



2 X 10-"gcm" 

Peri 



(10) 



where B, Q, Cs, p and pcri denote the magnetic field strength, 
the angular rotation speed, the sound speed, the central 
density of the cloud and critical density, respectively. If 
we Eissume the central density, critical density and tem- 

"gcm" 

and T = lOK of 
the molecular clouds, the above equation can be rewritten 
as follows: 



perature to be p = 2 x 10 ^Ogcm ^{n = 5 x 10^ cm ^) 
p^ri = 2 x 10"^^gcm-^(n = 5 x 10^° cm" 



n > 2.0 X 10" 



(img) 



1/2 



+ 0.25 



(11) 



Table shows the angular rotation speed (upper 8 clouds) 
and magnetic fiux density (lower 5 clouds) of severa l 
molecular clouds observ ed bv|Arquilla fc Go ldsmithI il98£tl , 
iGoodman et al. I Jl990h and lCrutcheJ il999l) as well as the 
magnetic fiux density and rotation speed necessary for frag- 
mentation expected from equation 1)11^ . For example, if 
the molecular cloud L1253 with angular velocity of Q — 
4.0 X 10~^'*s-'^ has a magnetic field weaker than B <3.1 /i 
G, the cloud fragments. On the other hand, if the molec- 
ular cloud L134 with magnetic flux density of _B = 11/iG 
has a rotation speed faster than f2 > 7.1 x 10"^* , the 
cloud fragments in the adiabatic accretion phase and has a 
possibility of forming a binary system. 



4.4 Effects of the Ambipolar diffusion 

In this subsection we discuss the effects of the ambipolar 
diffusion, which is not taken into account in this study. As 



discussed in paper I, the ambipolar diffusion has no appre- 
ciable effects as far as the isothermal collapse phase is con- 
cerned. However, the ambipolar diffusion plays an impor- 
tant role for high-density gases of n > 10^^"^^ cm"^ (see, 
e.g., iNakano et al. ll2002l . and the references therin). The 
ambipolar diffusion may change the outflows and fragmen- 
tation seen in our models seriously. 

First we examine the effects on the two types of out- 
flows, i.e., outer and inner outflows, shown in §3.2 and 3.3. 
The outer outflow emerges from the ring or bar structure, 
while the inner one does from the dense adiabatic core. The 
effect should not be serious for the outer outflow, since the 
density does not reach ~ 10^^ cm"'^ even at the base of the 
outflow. On the other hand, the inner outflow should weaken 
substantially since the magnetic held is coupled only weakly 
with the gas at the base. Thus the outer outflow will domi- 
nate if we include the ambipolar diffusion. 

Next we consider how the ambipolar diffusion effects 
the fragmentation. If the adiabatic core has an appreciable 
amount of angular momentum, it will fragment as shown in 
§4.1 and iDurisen et al. I il986h . On the contrary, the adi- 
abatic core will not fragment either if the core looses its 
angular momentum through the outflow and the magnetic 
braking, or if it has a small angular momentum at the mo- 
ment of formation. The ambipolar diffusion reduces the effi- 
ciency of the outflow and magnetic braking. Thus it renders 
to promotion of fragmentation. 

If the effect of ambipolar diffusion is restricted to the 
reduction in the angular momentum transfer, our conclusion 
that a cloud of a large Q/B fragments, is secure. This is be- 
cause the angular momentum transfer is already less efficient 
in our models showing fragmentation. The reduction in the 
angular momentum transfer will not change the result. On 
the other hand, the adiabatic core has only a small amount 
of the angular momentum at the moment of formation in 
the models showing no fragmentation. If the angular mo- 
mentum is small, the fragments will not survive but merge 
even if once they are formed. 



ACKNOWLEDGMENTS 

We have greatly benefited from discussion with Prof. M. 
Y. Fujimoto, Prof. A. Habe and Dr. K. Saigo. Numeri- 
cal calculations were carried out with a Fujitsu VPP5000 at 
the Astronomical Data Analysis Center, the National As- 
tronomical Observatory of Japan. This work was supported 
partially by the Grants-in-Aid from MEXT (15340062, 
14540233 [KT], 16740115 [TM]). 



REFERENCES 

Abt H. A., 1983, ARA & A, 21, 343 
Arquilla R., Goldsmith P. F. 1986, ApJ, 303, 356 
Basu S., Mouschovias T.C., 1994, ApJ, 432, 720 
Bate M., 1998, ApJ, 508, L95 

Bate M. R., BonneU I. A., Bromm V., 2002, MNRAS, 336, 
705 

Belloche A., Ande P., Despois D., Blinder S. 2002, A&A, 
393, 927 



12 M. N. Machida, T. Matsumoto, T. Hanawa and K. Tomisaka 



Binney J., Tremaine S. 1987, Galactic Dynamics (Princeton 

Univ. Pr., Princeton) 
Bodenheimer P., Burkert A., Klein R. I., Boss A. P. 2000, in 

Protostars and Planets IV (Tucson: University of Arizona 

Press) eds Marmings, V., Boss A.P., Russell S. S., p. 675 
Bonnell I., 1994, MNRAS, 269, 837 
Bonnell I., Bate M. R., 1994, MNRAS, 269, L45 
Bonnell I., Bate M. R., 1994, MNRAS, 271, 999 
Boss A. P., 1980, ApJ, 237, 866 
Boss A. P., 2002, ApJ, 568, 743 
Boss A. P., 2004, MNRAS, 350, L57 
Cha S.-H., Whitworth A.P., 2003, MNRAS, 340, 91 
Clarke C. J., Pringle J. E. 1991, MNRAS, 249, 588 
Cohen M., Kuhi LV., 1979, ApJL. 41,743 
Crutcher R. M. 1999, ApJ, 520, 706 
Duquennoy A., Mayor M., 1991, A&A, 248, 485 
Durisen R.H., Gingold R.A., Tohline J.E., Boss A.P., 1986, 

ApJ, 305, 281 

Dyck HM, Simon T, Zuckerman B., 1982, ApJL, 255, L103 
Goodman A. A., Batsien P. J., Puller G. A., Myers P. C., 

1990, ApJ, 359, 363 
Hosking J. G., & Whitworth A. P., 2004 MNRAS, 347, 3 
Joy A. H van Biesbroeck G., 1944. Publ.Astron.Soc.Pac. 

56,123 

Kroupa P., Burkert A., 2001, ApJ, 555, 945 
Larson R. B., 1969, MNRAS, 145, 271 
Machida M. N., Tomisaka K., Matsumoto T., 2004, MN- 
RAS, 348, LI (MTM04) 
Machida M. N., Matsumoto T., Tomisaka K., Hanawa T., 

2004, submitted to MNRAS (Paper 1) 
Masunaga H., Inutsuka S., 2000, ApJ, 531, 350 
Matsumoto T., Hanawa T., 2003, ApJ, 595, 913 
Miyama S. M., Hayashi C, Narita S., 1984, ApJ, 279, 621 
Nakamura P., Hanawa T., Nakano T. 1995, ApJ, 444, 770 
Nakamura F., Li Z.-Y. 2002, ApJ, 566, LlOl 
Nakamura P., Li Z.-Y. 2003, ApJ, 594, 363 
Nakano T., Nishi R., Umebayashi T. 2002, ApJ, 573, 199 
Norman M.L., Wilson J.R., 1978, ApJ, 224, 497 
Ohashi N., Hayashi M., Ho P.T.P, Momose M, Hirano N., 

1996, ApJ, 466, 957 
Richichi A., Leinert Ch., Jameson R., Zinnecker H., 1994, 

A&A, 287, 145 
Tamura M., Hough J. H., Hayashi S. S., 1995, ApJ, 448,346 
Tohline J.E., 1980, ApJ, 235, 866 

Tohline J. E. 1982, Fundamentals of Cosmic Physics, 8, 1 

Tomisaka K. 1996, PASJ, 48, 701 

Tomisaka K. 1998, ApJ, 502, L163 

Tomisaka K. 2000, ApJ, 528, L41 

Tomisaka K. 2002, ApJ, 575, 306 

Tsuribe T., Inutsuka S. 1999, ApJ, 523, L155 

Ward- Thompson D., Kirk J. M., Crutcher R. M., Greaves 

J. S., Holland W. S., Andre P., 2000, ApJ, 537, L135 
Whitworth A. P., Chapman S. J., Bhattal A. S., Disney M. 

J., Pongracic H., Turner J. A., 1995, MNRAS, 277, 727 



Collapse and Fragmentation of Molecular Clouds. II. 



Table 1. Model Parameters for Typical Models 



Model 


a 




A, 




1c,0 




-Bzc.O (>iG) 


f2c,0 (lO-* rad yr-l) 


M (Mq) 


L (10^ AU) 


M/MB,cri 


BL 


0.1 


0.01 


0.1 


0.2 


5 X 10^ cm- 


-3 


0.931 


1.26 


12.2 


6.82 


4.2 


CS 


0.01 


0.5 


0.1 


0.01 


5 X 10^ cm- 


-3 


0.295 


63.1 


20.6 


5.94 


UA 


DL 


1 


0.5 


0.1 


0.2 


5 X 10^ cm" 


-3 


2.95 


63.1 


28.7 


5.71 


1.4 
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a 



no 



£ob 



frag, mode survive 



1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 












0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.001 
0.005 
0.005 
0.005 
0.005 
0.005 
0.005 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.05 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 

1 
1 
1 
1 
1 
1 
1 
1 
1 

2 
2 
2 
2 
3 
3 
3 
3 
3 

1 (2) 
1 (2) 
1 (2) 
0.1(1) 


0.01 
0.01 
0.01 
0.1 
0.1 

1 



0.01 
0.02 
0.03 
0.04 
0.05 
0.1 
0.3 
0.5 
0.01 
0.02 
0.03 
0.04 
0.05 
0.1 
0.2 
0.3 
0.4 
0.01 
0.03 
0.05 
0.1 
0.3 
0.5 


0.01 
0.03 
0.05 
0.1 
0.2 
0.3 
0.4 
0.5 
0.1 


0.01 
0.03 
0.05 
0.1 
0.2 
0.3 
0.4 
0.5 


0.01 
0.03 
0.05 
0.1 
0.2 
0.3 
0.4 
0.5 
0.03 
0.4 
0.5 
0.6 


0.01 
0.03 
0.05 
0.1 
0.3 
0.4 
0.5 
0.6 
0.05 
0.1 
0.2 
0.3 
0.5 
0.3 
0.5 
0.5 



0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 



(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 
(0.2 



0.3 
0.3 
0.3 
0.001 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 



52 b 

52 
52 
52 
52 

c;2 



52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 
52 



52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

52 

5'- 

5- 

5- 

5- 

5- 
54 (5«) 
54 (56) 
54 (56) 
54 (56) 
54 (5«) 
54 (56) 
54 (56) 
54 (56) 



1.3 
2.9 
3.1 
3.5 
4.8 
5.8 
3.5 
5.1 
1.6 
1.9 
2.7 
3.7 
4.9 
6.0 
3.0 
3.4 
5.1 
2.1 
3.1 
4.7 
6.3 
3.5 
5.3 
2.9 
2.9 
3.7 
5.3 
5.8 
3.1 
3.4 
4.2 
5.1 
5.1 
5.3 
5.3 
5.7 
5.1 
4.5 
3.5 
3.5 
4.0 
5.2 
3.8 
3.7 
3.8 
3.7 
3.7 
4.0 
3.9 
4.8 
4.8 
5.1 
5.5 
4.3 
7.7 
6.7 
6.0 
6.4 
1.4 
5.9 
3.9 
4.4 
5.8 
4.2 
1.7 
3.4 
5.2 
4.8 
4.3 
4.5 
3.9 
3.1 



(1.4) 


3.9- 


3 


(6.4- 


2 N 







(0) 





17 


(0.16) 


(2.9) 


4.0- 
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The column Sob > Ear , B* , and Q* are the oblateness, the axis ratio, the square root of the magnetic pressure 
normalized by the thermal pressure [Bzc/iSiTC^pc)^^'^] and the angular speed normalized by the froefall 
timcscalc [flc/ {4:TtGpc)^^^] at core formation epoch of nc = 5 X lO-*^'^ cm-'^, respectively The column 'frag, 
mode' represents the fragmentation mode, and 'ring', 'bar' and ' — ' represent 'ring fragmentation model', 
'bar fragmentation model' and 'no fragmentation model', respectively. The column headed by 'survive' shows 
whether the fragments survive or merge after fragmentation as Y; survive to form a binary system, N; merge 
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Table 3. The Magnetic Flux Density and The Rotation Speed 
Necessary for Fragmentation 



Object 


n (10-" s-i) 


B (mG) 


L1174 


2.8 


< 1.3 


B35A 


3.4 


< 2.1 


L1253 


4.0 


< 3.1 


L1257 


5.6 


< 6.5 


B163 


6.6 


< 9.3 


TMC-2A 


7.1 


< 10.9 


B163SW 


9.7 


< 21.6 


L1495NW 


12.6 


< 36.7 


Taul6 


> 5.8 
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L134 


> 7.1 


11 


TauG 


> 8.5 


16 


W22B 


> 9.0 
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W49B 


> 9.7 
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Figure 1. The density (false color and contours) and veloc- 
ity distributions (arrows) for Model BL [(a, ui, j4^)=(0.1, 0.01, 
0.2)] are plotted on the z = plane at tc = 46 yr (a) , 182 
yr (b), 217 yr (c), and 227 yr (d) after core formation epoch of 



: 5 X IQio ( 



The time elapsed from the beginning of the 



calculation, the maximum density and the scale of the velocity 
vectors arc displayed in the upper area of each panel. The time 
after the adiabatic core formation is shown in the lower region of 
each ijaiicl. Levels of grid are shown in the upper left corner of 
the boundaries of subgrids. The thick line denotes the adiabatic 
core, namely, it means the contour of n = 5 x 10^" cm~^. 
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Figure 2. In each upper panel, the density (false color and 
contours) and velocity distributions (arrows) on x = plane are 
plotted for model BL [(«, cj, A^)={0.1, 0.01, 0.2)] at the same 
epochs of Fig. Thick lines denote the adiabatic core, namely, 
they mean the contour of n = 5 X 10^'^ cm~^. In each lower panel, 
t)z (false color and contours) and velocity vectors (vy , Vz : arrows) 
on the x=0 plane are plotted. The maximum speed in the z- 
direction (vz,max) is also shown in the upper part of the lower 
panels. Thick lines denote the boundary of Vz = 0. Inside the 
thick line, gas has an upward velocity [v^, > 0). 
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Figure 3. The magnetic field lines and outflow region are plotted 
for model BL [{a, lu, A^)=(0.1, 0.01, 0.2)] at the same epoch of 
Fig.lTlfd). Magnetic field lines (stream lines), adiabatic core (iso- 
density surface) and velocity vectors (arrows) on the z=0 plane 
are plotted in the left panel. In the right panel, the outflow re- 
gions (isovelocity surface of = 0) are also plotted beside the 
magnetic field lines, the adiabatic core and the velocity vectors. 
The box size is 40 AU. 
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Figure 4. The density (false color and contours) and velocity 
distributions (arrows) for Model CS [{a, w, A^)=(0.01, 0.5, 0.01)] 
are plotted on the 2 = plane at tc = 473 yr (a) , 879 yr (b), 
1344 yr (c), and 2134 yr (d). The contours, arrows, and notation 
have the same meaning as in Fig. [Tl 
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Figure 5. In each upper panel, the density (false color and 
contours) and velocity distributions (arrows) on y = plane are 
plotted for CS [(o, u), A^) = (0.01, 0.5, 0.01)] at the same epochs 
of Fig. 1^ In each lower panel, (false color and contours) and 
velocity vectors {vy, v^: arrows) on the x=0 plane are plotted. 
The contours, arrows, and notations have the same meaning as 
in Fig.[T] 
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Figure 6. The magnetic field lines and outflow region are plotted 
for model CS [{a, ui, A^) = (0.01, 0.5, 0.01)] at the same epoch 
of Fig. |4]{d). Magnetic field lines (stream lines), adiabatic core 
(isodensity surface) and velocity vectors (arrows) on z=0 plane 
are plotted in the left panel. In the right panel, the outflow regions 
(isovelocity surface of v^. = 0) are also plotted beside the magnetic 
field lines, the adiabatic core and the velocity vectors. The box 
size is 180 AU. 
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Figure 7. The density (false color and contours) and velocity 
distributions (arrows) for Model DL [{a,u>,A^) = (1, 0.5, 0.2)] 
are plotted on the ^ = plane at tc = 773 yr (a) , 1033 yr (b), 
1133 yr (c), and 1235 yr (d). The contours, arrows, and notations 
have the same meaning as in Fig. [Tl 
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Figure 8. In each upper panel, the density (false color and 
contours) and velocity distributions (arrows) on x = plane are 
plotted for Model DL [{a,uj,A^) = (1, 0.5, 0.2)] at the same 
epochs of Fig.|7| In each lower panel, v^, (false color and contours) 
and velocity vectors {vy, v^: arrows) on the x=0 plane are plotted. 
The contours, arrows, and notations have the same meaning as 
in Fig.Q] 
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Outer ojlllo^ 




Figure 9. The magnetic field lines and outflow region are plotted 
for DL [{a,u),A^) = (1, 0.5, 0.2)] at the same epoch of Fig. □ 
(d). Magnetic field lines (stream lines), adiabatic core (isodensity 
surface) and velocity vectors (arrows) on z=0 plane are plotted in 
the left panel. In the right panel, the outfiow regions (isovelocity 
surface of = 0) are also plotted beside the magnetic field lines, 
the adiabatic core and the velocity vectors. The box size is 300 
AU. 
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Figure 10. Pattern of fragmentation is shown on the magnetic 
flux density - rotation speed plane. The lower axis indicates the 
square root of the magnetic pressure (Bzc/VStt) normalized by 
the square root of the thermal pressure (y/^Pc) and the left 
ordinate does the angular speed (He) normalized by the freefall 
timescale (V^TrGpc), respectively. The upper axis indicates the 
parameter, a. Large symbols mean the fragmentation types: ring 
□ , bar A, both Q ^nd no fragmentation X are plotted against 
the initial values of B^c/iSnclpc)'^/'^ and nc/(47rGpc)^''^. SmaU 
symbols indicate whether fragmentation occurs o (open: ring frag- 
mentation, filled; bar fragmentation) or not +. These symbols are 
plotted against Bzc/{8-Kc^pc)^^'^ and nc/(47rGpc)^/^ at the end 
of the isothermal phase (nc = 5x10^" cm~^). Blue, red, and violet 
mean the regions of models of the ring, bar, and both fragmenta- 
tions for the initial parameter space for nc,o = 5 x 10^ cm~^. The 
two broken lines indicate the fragmentation region for rtc.o = 
5 X 10^ and 5 X 10"'' cm^'^, respectively. The dotted line indi- 
cates Qc/ {'iT^Gpc)^^'^ = 0.2. Four arrows represent the evolutional 
paths from the initial state (nc.o = 5 X 10^ cm""^) for models AS 
[{a,uj,A^) =(0.01, 0.01, 0.01)], BL (0.1, 0.01, 0.2), OS (0.01, 0.5, 
0.01) and DL (1, 0.5, 0.2), respectively. Thick line represents a 
curve of the magnetic flux - spin relation found in Paper I as 
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Figure 11. The angular rotation speeds normalized by the free- 
fall timescale are plotted against the oblateness (upper panel) 
and the axis ratio Ear (lower panel). These are values at the be- 
ginning of the adiabatic phase (rtc = 5 X 10^*^ cm^'''). The symbol, 
o, o, and + represent the ring, bar and no fragmentation, respec- 
tively. 
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Figure 12. The total (jtot), orbital (jorb) «m<i spin (jspn) specific 
angular momenta are plotted against the time after fragmentation 

for 6 models: (a) {a, a;, A^)= (0.01, 0.5, 0.01), (b) (0.1, 0.6, lO-^), 
(c) (1.0, 0.6, 10-'^), (d) (0.01, 0.6, 0.2), (c) (0.1, 0.5, 0.3), and (f) 
(1.0, 0.5, 0.2). Distance between each fragment (broken line) and 
fragments mass (dot-dashed line) is also plotted in each panel and 
the axes are displayed on the right-hand side of this figure. The 
upper three panels indicate the ring fragmentation models, while 
the lower three are the bar fragmentation models. The contour 
plot in each panel means the density distribution at about 1200 
years after the fragmentation. 



